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^ ' Abstract The Adler sum rule for deep inelastic neutrino scattering measures the isospin 
' of the nucleon and is hence exact. By contrast, the corresponding Gottfried sum rule for 
charged lepton scattering was based merely on a valence picture and is modified both 
by perturbative and by non-perturbative effects. Noting that the known perturbative 
corrections to two-loop order are suppressed by a factor l/N^, relative to those for higher 
moments, we propose that this suppression persists at higher orders and also applies to 
' higher-twist effects. Moreover, we propose that the differences between the corresponding 
radiative corrections to higher non-singlet moments in charged-lepton and neutrino deep 
^ ! inelastic scattering are suppressed by in all orders of perturbation theory. For 

^ I the first moment, in the Gottfried sum rule, the substantial discrepancy between the 
measured value and the valence-model expectation may be attributed to an intrinsic isospin 
asymmetry in the nucleon sea, as is indeed the case in a chiral-soliton model, where the 
discrepancy persists in the limit Nc oo. 
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1 Introduction 



Alone among the various sum rules of deep inelastic scattering (DIS) the isospin Adler 
sum rule P] has the special feature that its quark-parton model expression 



1 dx 



X 

1 



x,Q-')-Fr{x,Q') 
dx (u{x) — d{x) — u{x) + d{x)j 



4/, 



coincides with its QCD extension and receives neither perturbative nor non-perturbative 
corrections (for a discussion, see Ref. [2]). Moreover, this sum rule is supported by the 
existing neutrino-nucleon DIS data, which show no significant variation in the range 
2 GeV^ < Q2 < 30 GeV^ and give ^ 



rcxp 
-'A 



2.02 ±0.40 



(2) 



Though the error-bars are quite large, the precision could in principle be improved by 
future uN DIS experiments at neutrino factories (for discussion of such a program, see 
Ref.gl). 

Within the quark-parton model, the corresponding isospin sum rule in the case of 
charged-lepton-nucleon DIS has the form 



HQ') 
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If the nucleon sea were flavour symmetric, with u{x) = d{x), we would obtain the original 
Gottfried sum rule P|, Jg = |, in strong disagreement with the most detailed analysis of 
muon-nucleon DIS data, by the NMC collaboration, which gave the following result 



/g(Q2 = 4 Qgv^) = 0.235 ± 0.026 . 



(4) 



In contrast to the Adler sum rule, the original quark-parton model expression for the 
Gottfried sum rule is modified by perturbative QCD contributions, analyzed numerically 
at the a^-level in Ref.[Zj. These corrections turn out to be small and cannot be responsible 
for the significant discrepancy between Iq and the naive expectation of |. This discrepancy 
can be associated with the existence of non-perturbative effects in the nucleon sea, which 
generate light-quark flavour asymmetry, and lead to the inequality m(x, Q^) < d{x, Q"^) 
over significant ranges of the Bjorken variable x (for reviews, see Refs. jHl El 

In this paper we examine the QCD corrections to the moments of parton-model den- 
sities, for non-singlet neutrino and charged-lepton DIS, with the N = 1 moments cor- 
responding to the Adler and Gottfried sum rules, and comment upon a striking feature 
which they exhibit in the large- A'^c limit ^I] at the two-loop level. 
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2 Radiative corrections at large Nc 

First we present an analytical result for the two-loop radiative correction that was evalu- 
ated numerically in Ref.]?] and then comment on its structure as Nc —>■ oo. 



2.1 Analytical two- loop correction to the Gottfried sum rule 

Following Ref. [7j, we write the radiative corrections to the = 1 non-singlet charged- 
lepton moment of Eq. Q, in the case of light-quark flavour symmetry, as 

Ig = A{as)C^'\as) , (5) 
with an anomalous-dimension term 

AM - l + (6) 




where /?o and /3i are the first two scheme-independent coefficients of the QCD /5-function, 
namely 

/3o = (^Ca-Inf) (7) 



3 3 



with Np active flavours and Casimir operators Cp = {N^ — l)/{2Nc) and Ca = Nc, in the 
fundamental and adjoint representations of SU(Ai"c). 

The one-loop anomalous dimension vanishes and the leading correction from the two- 
loop result of Ref. |12j , confirmed in Ref. , has the form 

7f =1 = -4(C| - CfCa/2) [13 + 8C(3) - 12C(2)] ^ 2.557552 , (9) 

with two conspicuous features: 

• the appearance of ({2) = tt^/Q, which is absent from even non-singlet moments of 
the charged- lepton-nucleon structure function F2, and from odd moments of the 
corresponding neutrino-nucleon structure function, but occurs at odd moments for 
charged-lepton scattering, and at even moments for neutrino scattering, by analytic 
continuation in of results from QCD Feynman diagrams [TH] : 

• the distinctive non-planar colour-factor, {Cp — CpCa/2) = 0{N^), which exhibits 
an 0{1/N^) suppression at large- iVc, in comparison with the individual weights Cp 
and CpCai which are associated with planar two- loop diagrams that do not show 
this large- iVc cancellation at two loops ^21 for moments N > 1. Nor is there any 
sign of such large- iVc cancellation in the three-loop results of [Hj , obtained for even 
moments. 
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The second factor in Eq. (0) comes from radiative corrections to the coefficient function, 
of the form 



{£)N=1 



(10) 



with a vanishing one-loop term, C| 



{£)N=1 



ficient C' 



{e)N=i 



The scheme-independent two-loop coef- 



can be defined through the general non-singlet Mellin moment of charged- 



lepton-nucleon {£) DIS scattering 



dx 



c'(2).(+)(x,i)+c(2).{-)(x,iy 



„N-1 



(ii: 



taken at = 1, where the expressions for the functions C(2)'(-)(x, 1) and C(2).(+)(x, 1) 
were given in Ref.[T^ and confirmed later with the help of another technique in Ref.|17j. 
The "1" in the argument of these functions denotes the choice of renormalization scale 
yU^ = Q^, where /x^ is associated to the MS-scheme and the coupling as is evaluated at Q^. 

Explicit numerical integration of the = 1 moment of Eq. ()11|) gave the result [7] 



3.695C^ 



1M7CfCa 



(12) 



with a contribution from the colour factor CpNp which was consistent with zero, to the 
accuracy of that numerical work. At the time, the approximate emergence in Eq. fll2|) of 
the same non-planar structure {Cp — CfCa/^), already observed in the two-loop A^ = 1 
anomalous dimension coefficient of Eq. Q, went unremarked. Now we are able to derive 
an exact result, by comparing the charged-lepton moments pi|) with the corresponding 
non-singlet moments of the F2 structure function for neutrino-nucleon (z/) DIS, which can 
also be expressed in terms of the functions C^'^^'^~\x, 1) and C^'^^'^~^\x, 1), but now in the 
combination 

1 



2 Jo 



dx 



C(2).(+)(x,l)-C(^)'(-)(x, 1) 



X 



(13) 



To obtain an analytic expression for the correction (jlTj) to the Gottfried sum rule we 
remark that the A^ = 1 case of the moment ()13|1 corresponds to the Adler sum rule, which 



is free of radiative corrections. Hence, C- 



{u)N=l 



and by elimination of 



dxC(')'W(x,l) = / dxC(')'(-)(x,l) 



(14) 



we obtain 



2x3/ dxC^^^'^-\x,l) 



(15) 



Noting that the C^'^^'^~\x, 1) term in Ref . [T6] is explicitly proportional to Cf{Cf — Ca/2), 
we are left with a single integration over x, multiplied by this distinctive non-planar colour 
structure. Unlike the contributions from C^'^^'^^\x, 1), this integral is free of singularities 
as x — i> 1, and hence requires no regularization. The integrand involves trilogarithms, but 
elementary integration by parts reduces it to a regular integral whose integrand involves 
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nothing more complicated than the product of dilogarithms and logarithms. Maple then 
provided a speedy evaluation of the numerical coefficient of CpiCp—C aI'^ to 20 significant 
figures, for which we readily found a simple fit with a rational linear combination of 
the expected structures {1, C(2), C(3), C(4)}. Increasing the accuracy of integration to 30 
significant figures we confirmed, with overwhelming confidence, the analytical form 



{l)N=\ 



141 21 45 
-^ + ^C(2)-^C(3) + 12C(4) 



Cp{Cp - CaI2) (16) 
3.69439249494141137892516966638 Cf{Cf - Ca/2) , 



which validates the first 3 significant figures of the approximate terms of Eq. (fT2j) . obtained 
in Ref . [Zj by the far more difficult procedure of evaluating an integral in Eq. (fTT| that has 
three apparently distinct colour factors and requires delicate regularization at the singular 
endpoint, x = 1, of the C^'^'>'^^\x^ 1) function, interpreted as a distribution. 

We now interpret the vanishing of the one-loop corrections to the anomalous dimension 
and coefficient function of the = 1 non-singlet moment of charged-lepton-nucleon DIS 
structure functions as a simple consequence of the vanishing of all radiative corrections 
to the Adler sum rule and the absence of a non-planar one-loop diagram that distin- 
guishes charged-lepton scattering from neutrino scattering. As already remarked, this 
makes the two-loop anomalous dimension coefficient 7f^^^ and the two-loop correction 
^{pN-i gQ]^gj^g_j^^gpg^^g^|;_ The first place that scheme- dependence may appear is in 
the three-loop anomalous dimension coefficient 7^^^, which appears in Eq. (jH)) at order 
a^, albeit divided by [3q. This contribution is in the process of calculation (see for example 
Ref. [TH] ) . We expect its contribution to be small in the MS-scheme, for reasons discussed 
in Ref. [7j , based on experience of next-to- next-to-leading order fits jTHI to the data on xF^ 
in uN DIS from the CCFR collaboration. 

Moreover we offer our first conjecture, which is that the 6 possible colour structures in 
the three-loop term 7^^^ will occur only in those 3 combinations suppressed in the large- 

limit, namely Cl^Cp - Ca/2), CfCa{Cf - Ca/2) and Ci.(C^ - Ca/2)Nf. If this guess 
turns out to be wrong, then much of our subsequent discussion will become questionable. 
It should be noted that this conjecture applies exclusively to the A^ = 1 moment of the 
non-singlet charged-lepton structure function F2. We derive it from the wider hypothesis 
that the differences between non-singlet moments of F2 in charged-lepton scattering and 
neutrino scattering will continue to exhibit non-planar suppressions, beyond the two-loop 
order at which we have observed them. Then the suppression of 7^^^ in charged-lepton 
scattering at large Nc becomes a special consequence of the complete vanishing of radiative 
corrections to the Adler sum rule. 

We also note how quickly the two-loop corrections change their colour structure when 
one considers moments with A^ > 1. For example the ratio 



is forced to take the value -R^^^ = 1 at A^ = 1, by virtue of the vanishing of radiative 
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corrections to the Adler sum rule. But for = 2, we obtained from Ref.^T] the ratio 



R^=' = - . .r"":::::"' „ fi^ 



0.505931104 



^2 



5.4183241A'2 - AN^Np - 8.4480127 



which is negative and small in magnitude at large N^-, and also at A''^ = 3 with Np = 3 
active flavours, where it takes the value -R^^^ = —0.117197668. Moreover the magnitude of 
R2 continues to decrease very rapidly with the moment-number, N, because the integral in 
the numerator of Eq. p7|l has an integrand that is strongly suppressed as a; ^ 1. Similarly, 
we expect the currently known results for the charged-lepton anomalous dimension 7^, at 
several even values of A^, to give little guidance as to the eventual value at A^ = 1, which 
must be obtained by analytic continuation of a complete set of even-A^ results. 



2.2 Planar approximation, renormalons and corrections 

The limit ^ 00 and the l/Ag-expansion ^I] are known to be rather useful for analyzing 
the non-perturbative structure of QCD. Here we will use this framework to characterize 
our conjecture about the perturbative corrections and then seek a non-perturbative con- 
sequence. 

To do this, we use the planar approximation formulated in Ref . |20j . In this approx- 
imation one retains, after extracting an overall factor of Cp, only those terms at order 
{as/ir)" that contain the leading A^c behaviour for each possible power of Np. In the case 
of the order (a^/vr)"' contribution to the coefficient function of Eq. (fTUj) this prescription 
then amounts to selecting 

n-l 

Iplanar = E C^f =^ A^^ ^"^^ , (19) 

i=0 

where the Cn l^~^ are pure numbers. By definition, the planar approximation differs from 
reality by (at most) terms of order 1/A^. So far we have seen that C[^l^~^ = 0, since there 
is no one-loop correction to the coefficient function, and that C2^l^~^ = ^2^0^"^ = 0, since 
only the colour structure Cp{Cp — CaI'^) = ~\CpN~^ survives at two-loop order in this 
moment, because of the vanishing of all radiative corrections to the Adler sum rule and the 
appearance of a non-planar factor in the difference between charged-lepton and neutrino 
structure functions at two loops. Now let us analyze the consequences of the rather strong 
conjecture that the planar approximation (fTIHl also vanishes at all orders n > 2. 

In general, when it is non-vanishing, a planar approximation provides us with infor- 
mation in two distinct limits, namely in the large- A"c limit and also in the large- A^p limit. 
The intriguing link that it provides between these limits is underwritten by the way the 
large-order behaviour of perturbation theory is built by renormalon singularities, as was 
considered in QCD in the pioneering work of Ref. j|21j and reviewed in detail in Ref. [22] • 
This leads one to expect that the asymptotic behaviour of terms in perturbation theory 
in nth order is of the form C„ ~ KP^n^nl (where Pq is the first coefficient of the QCD 
/3-function) and so it is natural to develop perturbative coefficients as an expansion in 
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powers of (3o. The planar approximation is indeed polynomial in /?o and hence can be 
rewritten as 

n-1 

cT='Un.r = c,Y. ciT^'Po-'-'K , (20) 

i=0 

where again the Cl^l^~^ are pure numbers. This expansion is closely related to the proce- 



dure of naive nonabelianization (NNA) or large-/3o approximation proposed in Refs. 
in which one replaces Np by {IIN^ — 3/5o)/2 (for recent apphcations see Refs. [23 1^). 
The expansion of Eq. (j2(jp in Nc/jSo can be regarded as involving different numbers of 
effective renormalon bubble chains involving powers of /?o t22|, inserted in planar dia- 
grams 120] • There is a related expansion in Nf/(3q which is obtained by replacing Nc by 
{Z(5o + 2Nf)/11 [2111211 



n-1 



C^f Iplanar = E ^if = Vo "^"^iV]. , (21) 



=0 



and here again the Cn\^~^ are pure numbers. This expansion, which has been termed 
the "dual NNA" , has no direct Feynman diagrammatic interpretation, but turns out to be 
rather useful in making estimates of perturbative corrections to various physical quantities 
(see for example Ref.^25|). 



We now consider how the planar approximation is related to renormalon singularities. 
Following the work of Parisi j2Sj one expects that there will be singularities in the Borel 
transforms of QCD observables. We stress that we are focusing here on a coefficient 
function, say C, and ignoring any anomalous dimension part, since the latter will not 
contain renormalon effects [55]. C will have a Borel representation 

C{a) = r dze-'/''B\C]{z) . (22) 
Jo 

Here a = a^/vr and B[C]{z) is the Borel transform. The work of Parisi implies that one 
expects branch point singularities in z along the real axis at positions z = ±Zn where 
Zn = An/jSo, n = 1,2,3,4, . . .. Those on the positive real axis are referred to as infrared 
renormalons (IRn), and those on the negative real axis as ultraviolet renormalons (UV„,). 
Near each of these singularities one expects the structure 



^ ^^^,+ 0(l±^ (23) 



where the sum is over the contributions of various operators, and the 5i exponents depend 
on their anomalous dimensions. The large-order asymptotic behaviour of the perturbation 
theory will be determined by the dominant renormalon singularity nearest the origin, and 
its corresponding operator with largest 5i. The analysis has been carried out for the Adler 
e"'"e~-annihilation function, and for moments of the DIS structure functions Fi, F2 and F3, 
in Ref.(nni- UVi gives the dominant contribution for the Adler e"'"e~-annihilation function, 
and contributes, together with IRi, to the moments of DIS structure functions. The same 
dimension-six operator gives the dominant contribution to UVi in all the cases considered. 



6 



In the planar approximation one finds the exponent 




(24) 



m 2/3o 2A 



and one obtains the asymptotic large-order behaviour for the coefficient function of the 
A^th non-singlet moment of F2 



Only the overall constant Kp^ depends on the moment taken; the remaining n-dependence 
is universal [301 • Notice that in fact the same n-dependence also applies to the moments 
of Fi and F3 jSni. 

Our present conjecture is that the non-singlet moments of F2 in charged-lepton DIS 
and in neutrino DIS have essentially the same planar approximation, as a consequence 
of some generalization of the Cutkosky rules that were investigated to two-loop order in 
Ref.jl6j. One obvious consequence is that Ki = for the Gottfried sum rule, since clearly 
there are no corrections to the Adler sum rule. For higher moments the Kj^ [N > 1) will 
be nonzero, but very simply related. At n = 2 loops, one sees from Eqs. (fTTj) and (fT^ that 
both the i and u non-singlet F2 moments are dominated by C^'^^'^^\x, 1), at large N^. If it 
remains true beyond two-loop order that only the (-I-) component receives a contribution 
from planar diagrams, then one would expect that QC^^^l^^^^^j. = C'n^''^ Ipianar "with the 
factor of 6 simply resulting from the normalization of the Adler and Gottfried sums rules 
in the most naive quark-parton model. Not only would we expect GC^'^^ — C^'^^ to be 
suppressed by a factor of but also to decrease rapidly with the moment number, A^, 

as is the case at two-loop order. 

So far we have considered only the leading UV renormalon contribution. One may 
anticipate that there is an equally important IRi contribution, but to compute the corre- 
sponding 6 one would need the anomalous dimensions of twist-four operators contributing 
to the operator product expansion (OPE) for the non-singlet moments of F2, which are 
not known explicitly. The expectation would, however, be that the corresponding constant 
Kj^ would vanish for A^ = 1, and for A^ > 1 should differ by a factor of 6 for the u and i 
DIS moments. 

Since the leading OPE corrections to the moments of DIS structure functions 

are connected with the leading IRi renormalon (for a review, see Ref. |22j). we thus ex- 
pect higher-twist contributions to the Gottfried sum rule to be suppressed by a factor of 
as/{7iNc) ~ 1/(A^^ log(Q^/A^)) as Nc — >■ 00, relative to comparable effects in the Bjorken 
sum rules [2111211, because in the Gottfried sum rule a renormalon chain starts to develop 
only in a non-planar three-loop diagram, while in the case of the Bjorken sum rules it 
starts to develop in a two-loop planar diagram. 




(25) 



In the large- A^c limit one finds the asymptotic behaviour. 




(26) 
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3 The nucleon sea at large 



The previous discussion leads us to believe that the naive quark-parton model expression 
for the Gottfried sum rule, namely /g = |, is not modified by perturbative effects, or by 
their resummations as renormalon chains generating higher-twist effects, in the large- iVc 
limit. But in the real world, at Nc = 3, the experimental data of the NMC collaboration 
(see Eq. (j3])) show a very significant discrepancy from the naive expectation of |. 

There are several ways out of this puzzle. One is to say that = 1/9 is not small 

enough for our considerations to be relevant. Another is to say that the suppression 
to two-loop order was an accident that will not be repeated at higher loops. To our 
minds, the most interesting response is to allow that 1/9 may be a small enough factor 
to take seriously, and that such a suppression of radiative corrections may persist beyond 
two loops and hence be reflected in a suppression of higher-twist corrections, associated 
with IR renormalons. Then that leaves the failure of the naive Gottfried sum rule to be 
explained by an intrinsically non-perturbative flavour asymmetry of the nucleon sea that 
is inaccessible to renormalon analysis but should still be apparent in the Nc ^ oo limit, 
to which we have appealed in our perturbative conjectures and their resummations. 

It was interesting to learn from the authors of Ref. [32] that this is indeed the distinctive 
feature of a chiral-soliton model based on the work of Ref. [33]. Briefly, their large- A^^c 
picture, at a very low normalization point, around 0.6 GeV, is as follows. Isosinglet 
unpolarized distribution functions are large, since they give rise to sum rules that are 
proportional to N^, isovector unpolarized distribution functions appear only at next-to- 
leading order in 1/Nc, with the Adler sum rule satisfied in the form 

^/a = 1 = £^ dx (uix) - d{x)) (27) 

where the integrand at x > corresponds to a "constituent" quark contribution and at 
X < to an antiquark contribution coming from u{x) — d{x) = — (ll{—x) — d{—x)^. The 
failure of the Gottfried sum rule at large Nc is attributed to the integral 

1(3/g - 1) = - y° dx {u{x) - d{x)) = 1^ dx (u{x) - d{x)) = OiN^) (28) 

which measures the flavour asymmetry of the nucleon sea at this very low normalization 
point. Values of Iq between 0.219 and 0.178 were obtained for a range of constituent 
quark masses between 350 and 420 MeV, in fair agreement with = 0.235 ± 0.026 at 
= 4 GeV^. Note, however, that the NMC data are at a substantially higher momentum 
scale than can be accessed directly by the chiral-soliton model. For that reason, the authors 
also compared their predictions for u{x) — d{x) with the parton distributions of Ref. [SI] , 
which were initialized at a comparably low scale. Here too, they claim fair agreement. 

There are, of course, several other approaches to the problem of estimating the light- 
quark flavour asymmetry of the nucleon sea, based on meson-cloud models, instanton 
models and other considerations (see the reviews of Refs.|8.'i and the recent work in 
Ref . [35] ) . We have highlighted the results of the chiral-soliton model because it is based 
on the large- Ac expansion, used throughout this work. 
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4 Conclusions 



Within the large-iVc expansion we have made the following conjectures, based on rather 
limited two-loop input: 

1. Within the framework of light-flavour symmetry, the radiative corrections to the Got- 
tfried sum rule are suppressed by a factor l/N^, relative to the typical expectation 
0{{Ncas/n)'') ~ l/(log((5V^^))" at n loops. We base this on the facts that they 
vanish at the one-loop level and are merely of order {as/ir)^ ~ 1 / {Nclog{Q^ / A"^))^ 
at n = 2 loops. 

2. We expect the unknown three- loop anomalous- dimension coefficient 7^^^ to be re- 
stricted to only 3 of 6 possible colour structures, namely Cp{Cf — Ca/2), CfCa{Cp — 
Ca/2) and CpiCp - Ca/2)Nf. 

3. We expect the ratio of the non-singlet moments, with > 1, for the charged- lepton- 
nucleon and neutrino-nucleon F2 structure functions, to maintain the naive ratio 6:1, 
at large Nc, within the framework of light-quark symmetric perturbative QCD, after 
one discounts quark-loop terms involving Npd^'^'^dabclNc, which will contribute to 
the neutrino-nucleon moments. We have exposed the behaviour QC^^^ /C^^^ = 
1 + 0{1/Nl) for all > 1 at 72 = 2 loops and expect it to persist at higher loop 
orders in the quenched approximation, Np 0. 

4. Moreover, even at finite iVc, we expect this ratio to tend to unity at high moment- 
number A^, as is the case at two loops. 

5. We expect higher-twist corrections, of order l/Q^, to follow the same patterns and 
hence to be negligible in the Gottfried sum rule at large Nc- 

6. In attempting to reconcile this large- iVc perturbative picture with the significant 
discrepancy between the measured value for the Gottfried sum rule and the naive 
expectation of |, we note with interest the low-energy picture of the nucleon as a 
chiral soliton in the large- iVc limit, which leads to an intrinsically non-perturbative 
fiavour asymmetry of the nucleon sea [22] • We believe that current phenomenological 
analyses which incorporate a flavour- asymmetric sea as non-perturbative input, as 
for example in the most recent parton distributions of Refs.|3Bl EH EH EH], capture 
the essence of this situation, in a manner that cannot be achieved by radiative 
corrections, or by their resummations in the form of higher-twist effects. 

Acknowledgements We are most grateful to W. van Neerven and to C. Weiss, for 
clarifying remarks. DJB and ALK thank the members of IPPP, Durham, UK, for hospi- 
tality during the process of these studies. ALK's work was done within the program of 
RFBR Grants N 02-01-00601, 03-02-17047 and 03-02-17177. 

Note added in proof Shortly after we submitted our paper, an impressive de- 
termination of three- loop non-singlet splitting functions appeared in Ref . |40j . Using 
that work, we are now able to determine the three-loop anomalous-dimension coefficient 
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7^=^ = -2/o dxP4^)+(x), with Pji^+ix) given by Eq. (4.9) of Ref-gO]. To evaluate it, 
we note that the corresponding integral of the splitting function Pj^f>~{x) of Eq. (4.10) of 
Ref. [iU] vanishes and hence that 7^^^ = 2 Jq dx[P^f^^ (x) — P^f^^ (x)] indeed has the colour 
structure that we anticipated. Performing the integral analytically, we obtained 



+ 192C(2)C(3) 
- 192C(2)C(3) 



+ Cj_ 



290 - 248C(2) + 656C(3) - 1488C(4) + 832C(5) 
Hi + £!2a2)-ia!!5a3) + lf?C(4)-448C(5) 



9 3 ' 9 
^ 161.713785- 2.429260 A^F 

by systematic reduction of integrals of harmonic polylogarithms to Euler sums j41J with 
weights up to 5. This result was checked, to 30 significant figures, by numerical integration 
of an integrand involving products of dilogarithms, obtained after integration by parts. 
Within the framework of light-flavour symmetry, it leads to radiative corrections 



3/g 



1 + 0.035521 as/n - 0.58382 a^/vr^ for A^^ = 3 
1 + 0.038363 as/n - 0.56479 a^/vr^ for Np = 4: 



that are even smaller than those estimated in Ref. [7j, since the anomalous dimension terms 
of order cancel about 30% of the order contribution from the coefficient function. 
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